Synchronization, quantum correlations and entanglement in oscillator networks 
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Synchronization is one of the paradigmatic phenomena in the study of complex systems. It has 
been explored theoretically and experimentally mostly to understand natural phenomena, but also 
in view of technological applications. Although several mechanisms and conditions for synchronous 
behavior in spatially extended systems and networks have been identified, the emergence of this 
phenomenon has been largely unexplored in quantum systems until very recently. Here we discuss 
synchronization in quantum networks of different harmonic oscillators relaxing towards a stationary 
state, being essential the form of dissipation. By local tuning of one of the oscillators, we establish the 
conditions for synchronous dynamics, in the whole network or in a motif. Beyond the classical regime 
we show that synchronization between (even unlinked) nodes witnesses the presence of quantum 
correlations and entanglement. Furthermore, synchronization and entanglement can be induced 
between two different oscillators if properly linked to a random network. 



^ , * Correspondence to roberta@ifisc.uib-csic.es 
pLn ' Synchronization is a paradigmatic and well studiedphenomenon in many biological, physical and social systems 
[l|-l3|, also proposed as a tool in view of applications [J], but almost unexplored in the quantum regime. Interesting 
exceptions [5|-|7[ deal with entrainment, i.e. synchronization as the response to an external driver. At the microscopic 
level, spontaneous synchronization was recently predicted in nanomechanical oscillators [8, 9*] and observed [To| in 
the classical regime, when quantum fiuctuations are neg lected. Indeed a first approximation to a great variety of 



quantum systems, such as electromagnetic modes |lll - ll3l |. trapped ions [ij, ll5| or nanoelectromechanical resonators 



qua i 

[16| , is given by a set of coupled quantum harmonic oscillators, susceptible to experience spontaneous synchronization. 
' Beyond physical systems, there is an increasing awareness that quantum phenomena might play an important role 
" in terms of efiiciency of biological processes. Several examples [17Hl9f have shown that we may have to consider 
^ 1 quantum dynamics to explain biological phenomena. 

O^l A first step to characterize quantum spontaneous synchronization, considering quantum fluctuations and corre- 
lations beyond the classical limit, has been taken in Ref.jl^ where synchronization between one pair of damped 
,—1 ■ quantum harmonic oscillators has been reported. Most of the classical literature deals with self-sustained phase oscil- 
J> I lators modeled by Kuramoto-type models, or with identical nonlinear oscillators studied through the master stability 
formalism [s'] . Ref . [20^ takes a different view focusing on synchronization during the relaxation dynamics of different 
linear oscillators driven out of equilibrium and exploring the key role of dissipation. Indeed, depending on the way 
in which damping occurs, a pair of oscillators with different frequencies can show synchronous evolution emerging 
after a transient, as well as robust (slowly decaying) non-classical correlations [20'|. When several dissipative quantum 
. oscillators coupled in a network are considered dissipation can act globally or locally (in a node) and, depending on the 
' correlation length in the bath with respect to the size of the system, a variety of surprising phenomena are observed. 
In this work we show how synchronization can actually be induced by local tuning of one (even newly attached) 
oscillator of a generic (regular or random) network. Synchronous behavior emerges in the whole network or in a part 
of it and witnesses robust quantum correlations and entanglement. Stemming both from the structure of the system 
and from the form of system-bath coupling we further show the possibility to tune the system to configurations in 
which nodes do not thermalize and relax into a synchronous and non-classical asymptotic state. 
. The form in which dissipation occurs in a spatially extended system has deep consequences. We stress that the 
importance of symmetries present in the system-bath coupling has been recognized in many contexts. In classical 
systems, this fundamental issue was already discussed in the seminal work of Lord Rayleigh analyzing damping effects 
on normal modes in vibrating systems (2l| . Indeed, the role of dissipation to reduce detrimental effects of vibrations 
is fundamental in many areas of mechanical, civil and aerospace engineering 22]. On the other hand, in the context 
of quantum systems, symmetries in the coupling between qubits and the environment allow for decoherence-free 
subspaces 23] , entangled states preparation [2^ . [25| and dissipative quantum computing (26l . Wf\ . 

In the following we show that the distribution and form of losses through the network amounts to synchronous 
dynamics in spite of the nodes diversity, and robust quantum correlations. Furthermore, steady entanglement can be 
generated between unlinked nodes by properly coupling them to a network. 



FIG. 1: (A) Network of oscillators (represented by the network nodes) dissipating into separate baths (SB, represented by the 
gray circles surrounding the nodes). Links, representing couplings, have different strengths (lines thickness) and nodes have 
different natural frequencies (corresponding to different colors as given in the color bar). (B) Network of oscillators dissipating 
into a common bath (CB). (C) Network of oscillators with dissipation restricted to one node, local bath (LB). 



We consider generic networks of N nonresonant, coupled quantum harmonic oscillators, given by the Hamiltonian 



where q-^ = (gi, g^r) is the vector of canonical position operators and p are momenta, satisfying [qj,pj] — i (we take 
h = 1 throughout the paper), and Hm,n — ^^m^mn + Amn(l — Smn) is the matrix containing the topological properties 
of the network (frequencies Wm and couplings Xmn)- The eigenmodes of the system Q result from diagonalization of 
this Hamiltonian through the transformation matrix 

Any realistic model needs to include also environment effects [28l - l30j and, depending on system and bath correlations 
lengths, different forms of dissipation can be envisaged for an extended network (see Methods). All units can dissipate 
into separate identical baths (SB), Fig[T]A_, as cavity optical modes [l^llO]- Otherwise, if the coherence length of 
the environment is larger than the size of the system (here given by the spatial extension occupied by the oscillators 
network in Eq.([T])), all the nodes feel a 'similar' dissipation. This common bath (CB), FiglTj3, is known to create 
decoherence-free subspaces [1^ |3l| and asymptotic entanglement js^ - lsE! ]. A third, limiting, case (LB) in which a 
specific oscillator d dissipates much faster than any other node (FigUP) is also considered (local bath). 

One of the key insights of our work comes from noting that the coupling of real oscillators to the bath (taken here 
to be equal, 7) differs from those of eigenmodes. The latter are found to be jKm (with = X^n-^nm: '^m^ = 1 
and ~ J'dm) meaning that, except the SB situation, the eigenmodes have different decay rates. Then for CB and 
LB only the least dissipative eigenmode will survive to thermalization, thus governing the motion of all oscillators 
overlapping with it. It is then useful identifying the less dissipating normal mode with smallest effective coupling, k^, 
and also the following one, k,,, such that < Further, in the case of an Ohmic bath with cut-off- frequency 

larger than the frequencies of the system, the normal-mode damping rates simplify to the form = 7 



Knowledge of the normal modes of a complex network and of their dissipation rates (or effective couplings) allows to 
fully characterize a large variety of phenomena. Indeed this is a simple but powerful approach, even if diagonalization 
of the problem needs to be performed numerically except in a few (highly symmetric) configurations. By the diago- 
nalization matrix J- and system-bath interaction Hamiltonian we obtain the conditions to have a dominating mode 
during a transient. This mode dissipating most slowly, < \Kj\ Vj 7^ cr, is found either for CB or LB. Even more 
important, it is possible to identify normal modes completely protected against dissipation, so that some network's 
nodes do not thermalize. Indeed, a normal mode a is protected against decoherence if — 0. For a pair of oscillators 
interacting with a CB, this condition is accomplished only in the trivial case of identical frequencies 33', ,35] but this 
is not the case when more than two nodes are considered. We find that asymptotic synchronized quantum states can 
then be observed even in random networks where all nodes have different natural frequencies. 

Full characterization of the quantum state evolution of the network comes from moments of all orders of the oscillator 
operators qj and pj [30{ . For SB, average positions (and momenta) are characterized by irregular oscillations before 
thermalization (Fig. [5]^). On the other hand, for dissipation in CB and after a transient, regular phase locked 
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FIG. 2: First order moments for initial conditions — (gi) = (ga) = 1.0, (52) = 0.0, and vanishing momenta in the case of 
an open chain of three oscillators with uji = 1.2 aj2, tiJa = 1.8 aj2, non-vanishing couplings A12 = A23 — 0.4 ojI, temperature 
T = 10 u)2 (Boltzmann constant taken to be unity), 7 — 0.07 cj|, bath cutoff 50 ijJ2, for SB (A) and for CB (B). Time is scaled 
in units of the inverse of the reference frequency cj2. 

oscillations can arise, as shown in Fig. [2j3. Synchronization between detuned nodes can be found during a rather 
long and slow relaxation, like in the case of just one pair [23|. Further, the oscillations can remain robust even 
asymptotically if the condition Kg- = is satisfied. 

Beyond the classical limit given by average positions and momenta, let us now consider the full quantum dynamics 
stemming from the evolution of higher moments (see SI). At the microscopic level, quantum fluctuations also oscillate in 
time (even for initial vacuum states for which first order moments vanish at any time). This collective periodic motion 
is associated to a slow energy decay and witnesses the presence of robust quantum correlations against decoherence 
[29} . Our approach points to a wide range of appealing possibilities in quantum networks. In the following we show 
how a whole random network (or a part of it) can be brought to a synchronized state retaining quantum correlations 
via local tuning of just one of the nodes, or how two external oscillators can be linked to a random network leading 
to their entanglement and locked oscillations. 

Collective synchronization by tuning one oscillator 

Let us consider an Erdos-Renyi random and dissipative network [3] of oscillators with different node frequencies, 
links and weights (Fig. [T|) focusing on the relaxation dynamics of energy and quantum correlations. The node 
dynamics is mostly incoherent and even if initializing the network in a non-classical state, quantum correlations 
generally disappear due to decoherence [29]. Independently on the form of the network, for dissipation in SB, all 
nodes thermalize on a time scale 7"-'^ (see Eq. ^ and Methods). As anticipated before, this is not the case in the 
presence of a dissipation acting not-uniformly within the network. 

Common dissipation hath 

In presence of CB, an arbitrary network of N nodes can reach a synchronized state before thermalization if there is 
a weakest effective coupling n^. As a matter of fact just by tuning one of the node frequencies, aj„, even maintaining 
fixed the rest of the network frequencies {uji^y\ and its topology (A^ couplings) it is possible to decrease the weakest 
coupling Kg.. This means that an extra oscillator of properly selected frequency {uji^y} (like a synchronizer) can be 
added to a random network, even if weakly coupled, and it will lead to a collective synchronization of the whole 
system at some frequency (ilo-), generally different from Fig. [3] displays the average global synchronization and 
quantum correlations established in the network. Synchronization arises after a transient across the whole network 
by tuning one of the frequencies w^, to a particular value Coy , while it is not present when moving a few percent away 
from this value. Equivalently one could have tuned one of the couplings Xy^i . In the following we consider separately 
the case in which k„ is significantly smaller than the other effective couplings and the case in which it vanishes. 

Conditions for global synchronization follow from small ratio between the damping rates of the two slowest normal 
modes R — n^j k.^^ — >■ 0. Interestingly, this is a necessary but not sufficient condition for collective synchronization. 
This is due to the fact that the presence of a slowly dissipating normal mode into the system needs to be accompanied 
by a significant overlap between this mode (Q^), or virtual oscillator, and all the real ones (qi, ...^qn). An analytical 
estimation of the synchronization time is found taking into account both the importance (overlap with individual 
oscillator) and decay of few normal modes in the system. The time for oscillator j to start oscillating at the less 
damped frequency fla- reads: 



(2) 
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FIG. 3: (A) Time evolution of synchronization S, and (B) quantum correlations quantified by the discord {S) x 10^ , when 
varying one node frequency Results are shown for a random network (connection probability p — 0.6) of 10 oscillators. 
Frequencies of nodes are sampled from a uniform distribution from O.Qojq to 1.2ajo and couplings from a Gaussian distribution 
around — O.Iojq with standard deviation O.OSojq. Environment parameters (here and in the following figures) are 7 = O.OltJo, 
T = lOuJo and cut-off frequency A = 50a;o (see Methods). Time units are chosen so that too = 1. Collective synchronization S 
and (averaged and filtered) discord S (see Methods) are obtained considering all oscillator pairs of the network. Dashed line 
identifies the frequency uu^ for which Kc — 0. Continuous line in (A) corresponds to the estimated synchronization time tsync- 



maximizing over all modes k different from the slowest one cr (see SI). Collective synchronization time corresponds 
to tsync = i.e. when even the last oscillator joins the synchronous dynamics dominated by the less 

damped mode. Then phase-locking in the evolution of all oscillators, namely in their (all order) moments, can arise 
before thermalization, when there is significant separation between largest time decays r^,ro-, and overlap between 
slowest normal modes and each system node. Global network synchronization (see Methods) obtained from the full 
dynamical evolution and the estimated synchronization time tsync are in good agreement, as seen in Fig. For the 
same network, the ratio i?, the collective synchronization S and mean discord {5) at long times can be seen in Fig. 
SI (see SI). 

We now look at the quantumness of the state in presence of collective synchronization. Generally decoherence 
is independent on the specific features such as the oscillation frequency in a system [36]. Still, synchronization 
is a consequence of a reduced dissipation in some system mode and indeed witnesses the robustness of quantum 
correlations, as evident form the average discord {5) in the network represented in Fig|3j3. After a transient dynamics 
in which the couplings in the network create quantum correlations [16|, even when starting from separable states, 
discord does actually decay to small values for lj^ different from Co^ (non-synchronized network) while it maintains 
large values for the case of a properly tuned node (w^ ~ Wu)- 

The case cj„ — oJv, leading to Kq- = 0, needs special attention. After a transient all the nodes will oscillate at a locked 
common frequency, the one of the undamped normal mode which we call a frozen mode. As before (Eq.([5])), 
the possibility to synchronize the whole network also requires a second condition, namely that the undamped mode 
involves all the network nodes. (The case in which the latter condition applies only to some nodes is discussed below.) 
When both the conditions 

N 

/tg^ = ^ Tka- = , and Tka ^ Vfc (3) 

fc=i 

are met, there is a frozen normal mode linking all oscillators. This leads to collective synchronization in the whole 
network and allows for mutual information and quantum correlations remaining strong even asymptotically, being 
orders of magnitude larger than for the fully thermalized state, when synchronization is not present (Fig|3]). The 
undamped mode gives actually rise to a decoherence-free dynamics for the whole system of oscillators where quantum 
correlations and mutual information survive. 

The phenomena above are found for nodes dissipating at equal rates into a CB, while in the presence of N indepen- 
dent environments (SB) all oscillators thermalize incoherently, synchronization is not found, and decoherence times 
for all oscillators are of the same order. As a final observation we mention the special case in which the center of mass 
of the system is one of the normal modes; then there will be a large decoherence-free subspace (corresponding to the 
other — 1 modes) but no synchronization will appear for a CB. 



Local dissipation bath, 



Common and separate baths correspond to two extreme situations in which all oscillators have equivalent interac- 
tions with the environment (s). We now consider the case of a local bath, as a limit case in which one oscillator is 
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FIG. 4: (A) Synchronization factors Sc and (B) average discord Sc x 10^ evaluated for linear 3-node motifs (hence the 
subindex C) in a random network (connection probability p = 0.6) of 15 oscillators (shown in C). A tuned non-dissipative 
motif Ci composed by the three nodes (a — c — 6) is compared with another equivalent non-tuned motif C2 composed by nodes 
(d — / — e). Frequencies in the network are sampled from a uniform distribution from ujq to l.Sujo, and couplings with a Gaussian 
distribution around —OAluq with standard deviation O.OSojq- In order to avoid dissipation in the {a — c — b) motif we have set 
LJc ~ 1.51 ujo, being Xac = —0.09 luq and Xtc ~ —0.11 ljq. Time units such that ujo = 1. 



dissipating stronger, FigHp. A frozen normal mode cr must not overlap with the dissipative oscillator (labeled by d) 
while involving all the other nodes {J-ia 7^ V i 7^ d). Then, synchronization of the whole network (except for the 
dissipative oscillator) arises. This occurs when 

J^da = , with J^dj 7^ V (4) 

meaning that the undamped mode a involves a cluster of oscillators not including the lossy one. We find synchroniza- 
tion and robust quantum effects across the network as for CB, with the difference that for LB the dissipating node is 
now excluded (further details are discussed in SI). 



Synchronization of linear motifs 



The possibility to synchronize a whole network, in presence of different dissipation mechanisms, just by tuning one 
local parameter opens-up the perspective of control that can be explored considering the dynamical variation of a 
control-node frequency. In particular we find similar qualitative results both for random networks and for disordered 
lattices consisting of regular networks with inhomogeneous frequencies and couplings, being the latter largely studied 
in ultracold atomic gases (37j . Local tuning to collective synchronization is not only a general feature of different kind 
of networks but can also be established in motifs within the network. As we show in the following, the system can be 
tuned to a partial synchronization, involving some nodes of a network independently of the rest of it. Indeed, even if 
the whole system is coupled to a CB, we can identify the conditions for having a synchronized cluster, like the 3-node 
linear motif in Fig. 01 Two non-directly linked nodes a and b of the motif are then asymptotically synchronized 
through another one, here c, and this leads to a common oscillation dynamics along the whole motif a-c-b. The 
condition for synchronization of a cluster, namely its dependence on a frozen normal mode, reads 

-2 = -1 (5) 



with r^CT frequency of the frozen mode (see details in SI). 

This case is an example of the general result stating that given any network, a part of it (in our case a linear motif, 
Ci) can be synchronized by tuning one of its components, for instance a frequency or coupling of the motif. A key 
point is that this is independent of the frequencies and links of the rest of the network, provided the motif is properly 
embedded in the network. The links between Ci and the rest of the network should satisfy 

+ (7^^^) + Ac, - , Vj . (6) 



This is equivalent to saying that a synchronized motif with robust quantum correlations can preserve these features 
when linked to an arbitrary network, if some constraints on the reciprocal links are satisfied. For instance, each node 
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FIG. 5: (A) Energy evolution and (B) entanglement (logarithmic negativity) between two nodes with identical frequency 
ujo (we call these nodes a and b and are plotted in red in the network displayed in panel (C)). The network is a random one 
(connection probability p = 0.6 of 15 oscillators and same frequency and couplings distribution as in Fig. 3] We compare the 
situations in which the couplings from the red nodes to the rest of the network (they are directly connected to other nodes 
called c and d) are properly balanced in order to avoid dissipation (Aac = \bc = — O.lStJo and Xad = Xbd = — 0.12cJo) with 
the case when this balance is perturbed (Aac + 0.04a;o and Xad + 0.04ajo) .The third line in panel (B) shows the entanglement 
between other two arbitrary oscillators in the situation in which a and b are balanced. Time units such that too = 1 



of the synchronized motif needs to share with the rest of the network more than one link. In Fig. 2] we compare the 
behavior of two linear motifs of a large network, where a first motif Ci is synchronized, satisfying Eqs.([5l)-(l6]) while 
the second one C2 is not. After a transient a frozen mode tames the dynamics of Ci, which then shows a synchronous 
evolution and robust correlations. It can also be shown that quantum purity and energy reach higher values of a 
stationary non-thermal state. This is compared with the non-synchronized motif C2 whose dynamics quickly relaxes 
to a thermal state. The case of a three oscillator chain is an example showing the possibility to tune synchronization 
and quantum effects in a motif within the network when the proper link conditions are satisfied. 



Entangling two oscillators through a network 



The same technique discussed in the previous section can be applied to the case in which we aim to synchronize 
few, even if not directly connected, elements of a network. But this does not mean that any set of nodes can be 
synchronized asymptotically. In fact, we find that for a CB we can only synchronize two different and not directly 
linked {Xat = 0) oscillators if we synchronize with them also other intermediate linked elements (like in the linear 
motif example. Fig. 2]) or when these two oscillators are identical, like we discuss in the present section. 

We consider the case of two identical oscillators (i.e. with uja = i-Ob) prepared in a separable state, with some local 
squeezing. They are not directly coupled {Xat = 0) but are connected through an arbitrary network. In general they 
will dissipate reaching the thermal state, but with the proper conditions we find an important result: because their 
frequencies are identical it is possible to construct a frozen normal mode involving only these two nodes, given by 
Qa = Taa Qa + J^bcr <Z6, with Taa^^ha 7^ and this Can be obtained, for instance, by properly attaching them to the 
network. In other words, by tuning their couplings to the network it is possible to have both oscillators relaxing onto 
a frozen mode, so that they will be synchronized and will keep a higher energy than otherwise. Most importantly, in 
this case entanglement can actually be generated between oscillators initially in a separable state and remains high 
asymptotically. 

In order to entangle the oscillators, their coupling to the rest of the network needs to fulfill the condition 

^ Tka\k, = 0, (7) 
h—a^h 

(similar to Eq. (jH])), achieved by proper tuning of coupling strengths of the active links (j) with the rest of the network 
Aaj, Af,j. In Fig. [5]4. and B we show the evolution of energy and entanglement of the oscillators a and h when linked 
to a random network. As we see in Fig. [Sp there is not direct link between a and h (Xgh — 0) and the whole system 
dissipates in a common environment. The case where the oscillators a and h are coupled to the network following 
the prescription ([7]) is compared to another case in which their links are not properly balanced (we slightly change 
the coupling strengths) . Both energy and entanglement are shown to be sensitive to the structure of the reciprocal 
links and the possibility to actually bring the added nodes into an entangled state that will survive asymptotically is 
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guaranteed by Eq. ([7]). The importance of this result is twofold: in terms of applications it shows that it is possible 
to dynamically generate entanglement between two non-linked nodes embedded in a random network by tuning their 
connections to it, and on the other hand it enlarges the scenario for asymptotic entanglement generation through 
the environment. It is known that large entanglement can be generated between far oscillators during a transient 
due to a sudden-switch or through parametric driving (38i] . On the other hand, a common environment leads to 
entanglement between a pair of spins [23| or oscillators [33| . 

DISCUSSION 

Our results on synchronization in dissipative harmonic networks and its optimization give a flavor of all the possi- 
bilities that show up once the mechanism behind the phenomenon is understood. At difference from widely considered 
self-sustained non-linear oscillators, here we focus in a linear system showing how synchronization can emerge after 
a transient for dissipation processes introducing inhomogeneous decay rates among the systems normal modes. A 
synchronous oscillation is predicted, for the first time, either in a long transient during relaxation to the equilibrium 
state or in a stationary non-thermal state. We considered the most significant examples of correlation length of the 
bath larger than the system size (CB) and of a node of the network more strongly exposed to dissipation (LB), 
displaying synchronous dynamics. On the other hand, for independent environments (SB) on different nodes the 
resulting dynamics remains incoherent even when increasing the strength of the reciprocal couplings in the network. 

The presence of synchronization in the whole or a part of the network witnesses the survival of quantum correlations 
and entanglement between the involved nodes. This connection between a coherent oscillation in the network and 
its non-classical state is a powerful result in the context of complex quantum systems, considering the abundance of 
this phenomenon. Indeed, the condition underlying synchronization provides a strategy to protect a system subspace 
from decoherence. Our discussion and methodological approach are general, but we show specific consequences of our 
analysis, such as global or partial synchronization in a network through local tuning in one node (synchronizer) as 
well as the possibility of connecting two nodes (not linked between them) to a network and synchronize and entangle 
them, even starting from separable states. Even if the reported results refer to random networks, our analysis applies 
to generic ones, also including homogeneous and disordered lattices and do not require all-to-all connectivity. 

In some sense, tuning part of a network so that the rest of it reaches a synchronous, highly correlated state can be 
seen as a kind of reservoir engineering, where here the tuned part of the network would be a part of the reservoir. This 
is to be compared with recent proposals of dissipative engineering for quantum information, where special actions are 
performed to target a desired non-classical state [23 - [27| . In the context of quantum communications and considering 
recent results on quantum Internet [39,, [40| , our study can offer some insight in designing a network with coherent 
information transport properties. Furthermore implications of our approach can be explored in the context of efficient 
transport in biological systems [l^, [l^ . An interesting methodological connection is also with transport through 
(classical) random networks [4lj. On the other hand, our analysis, when restricted to the classical limit, also gives 
some insight about vibrations in an engineering context, providing the conditions for undamped normal modes and 
their effect 

METHODS 

Interaction with the environment 

The Hamiltonian of the system Hs as defined in Eq. ([T|) is diagonalized in the basis of its eigenmodes Q = ^"^q 
yielding — J-^T-LJ-. In a microscopic description with independent oscillators modeling the environment, the 
system-bath interaction Hamiltonian for SB takes the form 

N oo 

Hf"" = -7 E '^'"^^"^ ' ^itl^ ^^"^ = E ^"^i"^ (8) 

m— 1 a—1 

being 7 the system-bath coupling strength (explicitly shown for the ease of understanding) , the position operators 
for each environment oscillator a (representing for instance a vibrational mode, or an optical one, etc..) of the bath 
in which the network unit (to) is dissipating. As explained in the main text, this situation occurs when the 
coherence length of the environment is smaller than the spatial extension of the system. Thus each oscillator dissipates 
to its own heat bath. In the opposite case, a common bath is seen by all oscillators, resulting in an interaction 
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Hamiltonian 

N 



Hf" = -lYl 9™^, (9) 



m— 1 



and actually involving only the average position (here the center of mass) of the network, Fig. [Tj3. Notice that in the 
eigenmodes basis 



ijf ^ = ~l'^H,nQmB , with K„ = ^ Tnra- (10) 
m n 

The effective couplings Km are different and determined by characteristics of the network such as topology, coupling 
strengths, and frequencies, as encoded in the diagonalization matrix F. This is in stark contrast to the case of identical 
SB ([U where all normal modes have equal effective couplings to the baths (to see this, notice that we can transform 
the bath operators Xq™'' to a new basis which exactly cancels the transformation J^; these new 'oscillators' can be 
shown to have the same statistical properties as the others, thus resulting in equivalent heat baths). 
Finally, the case of a given node d dissipating much faster than any other is modeled by 

Hf"" = -IQdB. (11) 

This local bath (LB) situation does also lead to non- uniform environment interaction in some of the normal modes 
with effective couplings Km'- 

ijf ^ = -'^'^KmQmB , with K„i Jdm • (12) 

m 

All the mentioned situations of separate, common and local bath can be described by different master equations 
[28| for the evolution of the network state. 



Master equation 



A standard procedure allows to obtain from the total Hamiltonian of system, bath and reciprocal interaction the 
evolution of the reduced density matrix for the state of the system, in our case a network of different oscillators. After 
a (post-trace) rotating wave approximation, the master equations in the weak coupling limit for separate, common, 
and local baths are in the Lindblad form, guarantying a well-behaved system d yna mics. These equations (given in 
SI) are obtained by generalization of the problem of a pair of coupled oscillators [20| ■ For the purpose of our analysis 
it is interesting to consider the master equation in the normal-mode basis {[Qk, Pi] = iSki)- 

- ([Q"' /^w}] ~ [^"' w»' '°(*)}]) + 

71 

1 



+ ( [Q„, [g„,p(t)]] - -j[P„, [Pn,p{m ) (13) 



where r2„ are the normal-mode frequencies of Hs and the damping and diffusion coefficients, for an Ohmic bath with 
spectral density J(a;) = {2^ /tt)u]Q{K — lu), are, at temperature T far away from the frequency cutoff A: r„ — 
and Dn — k^j70„ coth(^), being k„ the effective couplings (see Eqs. pUj) and ((HI)). Boltzmann constant is taken 
to be unity. With the appropriate definition of the couplings this equation is valid both for common and local bath 
while for CB we have: r„ — 7 and Z)„ = 7r2„coth(^) i.e. we obtain the same damping coefficient for all normal 
modes. The main differences in the models of dissipation here proposed reside in these expressions for the master- 
equation coefficients that will produce different friction terms in the equations of motion determining collective or 
individuals behaviors (See SI). We stress that the choice of this master-equation representation is not critical for our 
main conclusions, as also shown in 12011. 
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Synchronization factor 

Synchronization between two time series f(t) and g{t) can be characterized by a commonly used indicator, namely 
Cf^g{t,At) = 5f6g/ \JYp~5^ where the bar stands for a time average f = Jt^^*^ dt' f{t') with time window At 

and Sf = f — f. For 'similar' and in phase (anti-phase) evolutions C ^ 1 (—1), while it tends to vanish otherwise. As 
a figure of merit for global synchronization in the whole network we look at the product (neglecting the sign) of the 
indicator for all pairs of oscillators in the system. When the time series correspond to positions second moments we 
have S = ni<j |C^^2^^^2^ |. This collective synchronization factor S can reach unit value only in presence of synchronous 
dynamics between all the pairs of oscillators in the network. 

Mutual information, discord and entanglement 

The total amount of correlations present in a bipartite system can be measured by the mutual information I = 
Sa+Sb — Sab with Si the entropy of the reduced system i = A,B and Sab the total entropy. It can be decomposed in a 
purely classical part Ja-b = min[5A - SA\Ef], with SA\Ef = J2iPiSiSA\E^): Pi = T^^AB{EPg) and Qaie^ = Efg/pi 
(where i labels the possible outcomes of a general measurement with operators acting on B occurring with 
probability pi and resulting in a density matrix gi^j^s), and a quantum part which is just the difference 5a:B = 
I — Ja:B^ the quantum discord (42| - |4^ . As a quantifier for entanglement we have used the logarithmic negativity 
-Ejv = max(0, — lnz^_), with the smallest symplectic eigenvalue of the partially transposed density matrix [i^. We 
quantify correlations between pairs of (linked or unlinked) nodes by applying these definitions to different pairs of 
oscillators in the network. Average quantities ((.)) for all pairs in the network can be considered to assess globally 
the network. 
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